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The complex elastic compliance and dielectric susceptibility of (Nao.5Bio.5)i-iBaa:Ti03 (NBT- 
BT) have been measured in the composition range between pure NBT and the morphotropic phase 
boundary included, < a:: < 0.08. The compliance of NBT presents sharp peaks at the rhombohe- 
dral/tetragonal and tetragonal/cubic transitions, allowing the determination of the tetragonal region 
of the phase diagram, up to now impossible due to the strong lattice disorder and small distortions 
and polarizations involved. In spite of ample evidence of disorder and structural heterogeneity, the 
R-T transition remains sharp up to a:: = 0.06, whereas the T-C transition merges into the diffuse and 
relaxor-like transition associated with broad maxima of the dielectric and elastic susceptibilities. An 
attempt is made at relating the different features in the anelastic and dielectric curves to different 
modes of octahedral rotations and polar cation shifts. The possibility is also considered that the 
cation displacements locally have monoclinic symmetry, as for PZT near the morphotropic phase 
boundary. 
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I. INTRODUCTION 

The phase diagrams of the solid solutions of the fer- 
roelectric perovskite Nao.sBio.sTiOa (NBT) with other 
tetrahedral perovskites are actively studied, mainly with 
the intent of obtaining compositions close to a mor- 
photropic phase boundary (MPB) between rhombohe- 
dral (R) and tetragonal (T) phases, with easily reori- 
entable ferroelectric domains and high electromechanical 
coupling. In fact, NBT based perovskites are among the 
possible lead-free piezoelectric materials that would sub- 
stitute the widely used PbZri_2,Tia;03, when new envi- 
ronmental legislations will limit the use of Pb containing 
materials. 

The phase diagram of (Nao.5Bio.5)i-2;Baj;Ti03 (here 
abbreviated as NBT-BT or BNBT-lOOx) is still uncer- 
tain in the interesting region a; < 0.1 including the MPB, 
particularly regarding the nature and extension of the 
presumably almost antiferroelectric (AFE) phase that 
becomes ferroelectric (FE) on further cooling. Also the 
temperature evolution of the structure of undoped NBT 
is the object of contradictory interpretations.^ The main 
difficulties in determining the structure of NBT by neu- 
tron or x-ray diffraction are due to the very low distor- 
tion in both the tetragonal and rhombohedral phases, 
with a cell distortion from cubic of < 0.2% and ~ 0.4% 
respectively,- and to the structural disorder in the Na/Bi 
sublattice, which is averaged out and compensated by 
anomalously large thermal factors in the analysis <^ In 
addition, x-ray diffraction is little sensitive to the rota- 
tions of the O octahedra involved in the various struc- 
tural transformations, due to the low atomic number of 
O. Cation disorder is enhanced in the solid solutions with 
other perovskites; it is responsible for the diffuse nature 
of the lower temperature transitions in the dielectric mea- 
surements, sometimes considered of the relaxor type, and 



possibly for the very broad temperature range of coexis- 
tence between T and R phases, with characteristics that 
vary from study to study. There is also scarce repro- 
ducibility of the transition temperatures and of the as- 
sociated anomalies in the dielectric and other physical 
properties, which generally exhibit very broad anoma- 
lies with large hysteresis between heating and cooling. 
These uncertainties may also be connected with differ- 
ences from sample to sample in the local cation ordering 
and non perfect stoichiometry in the Na/Bi sublattice, 
due to preferential loss of Bi during sintering. 

The elastic constants should provide clear information 
on the occurrence of transformations between C, T and R 
phases, which are of ferroelastic nature, but to our knowl- 
edge only few studies of the acoustic properties of NBT 
exist which have not an obvious connection with the 
known phase transitions. Here we present an extensive 
study of (NBT)i_2;(BT)^ with 0<x < 0.08 by combined 
anelastic and dielectric spectroscopies. 



II. EXPERIMENTAL 

Powders of the composition (Nao.5Bio.5)i-a;Baa;Ti03 
with X = and x = 2, 3, 4, 5, 6, 8 mol% were pre- 
pared following the mixed-oxide method. Stoichiomet- 
ric amounts of NajCOa (Merck 6392), BizOs (Aldrich 
223891), BaCOa (Merck 1714), TiOs (Degussa P 25) 
were mixed in ethanol with zirconia milling media for 
48 h, dried and sieved. For each composition the calcin- 
ing and sintering temperatures were systematically in- 
vestigated in order to improve the cold consolidation be- 
havior and the final density. After heat treatment at 
700 — 800 °C for 1 h and further milling and drying, the 
powders were isostatically pressed at 300 MPa and sin- 
tered at 1150 °C for 2 h. In order to avoid the loss of Na 
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and Bi, which is significant at temperatures over 1000 °C, 
most of the sintering processes were carried out with the 
samples placed on a Zr02 disc, covered with an AI2O3 
crucible and sealed with NBT pack. The x-ray diffraction 
analysis evidenced the pure perovskitic phase in all the 
samples investigated. The microstructure of the sintered 
samples was observed on polished and etched surfaces 
by means of scanning electron microscopy (Leica Cam- 
bridge Stereoscan 360). Depending on the composition 
(at the same sintering conditions) the microstructures re- 
sult more or less homogeneous, with some large, squared 
grains in a matrix of smaller grains; among those sintered 
in the NBT pack powder, pure NBT showed the largest 
grain size (5 — 20 ^m), which reduced upon addition of 
barium, the lowest mean size being found with x = 0.08 
(0.5 /^m). 

We prepared discs (diameter 22 mm, thickness 2 mm) 
for the dielectric measurements and blocks (60 x 7 x 
7 mm^) subsequently cut into bars (40 x 4 x 0.6 mm'^) for 
both anelastic and dielectric experiments. The faces of 
the discs were ground, silver electroded and poled with 
an electric field of 3 kV/mm in a silicon oil bath at 120 °C 
for 30 min, while the bars were tested unpoled. 

The dielectric measurements were carried out on poled 
and unpoled discs and bars in a frequency range between 
200 Hz and 200 kHz using a Hewlett-Packard LCR meter 
(HP 4284A) with a four wire probe and an ac driving 
signal level of 0.5 V/mm. The dielectric permittivity 
e = e' — ie" was obtained from the measured values of 
capacitance and loss tan J = e"/e'. The measurements 
were made on heating/cooling at 1 — 1.5 K/min between 
300 and 570 K in a Delta Design climatic chamber model 
9023A. 

The complex Young's modulus E = E' + iE" was mea- 
sured by electrostatically exciting the flexural vibrations 
of the bars, which were suspended on thin thermocouple 
wires in vacuum. It was possible to excite the 1st, 3rd 
and sometimes 5th flexural modes, whose resonating fre- 
quencies are in the ratios 1 : 5.4 : 13. The fundamental 
frequency is given by^ 




where I is the length, h the thickness, p the density the 
bar; for most samples it was / ^ 1.4 kHz. From the above 
formula it is possible to extract the absolute value of the 
Young's modulus. Such a value is affected by an error 
that may easily be of the order of 10%, due to uncertain- 
ties in the sample dimensions, porosity (the density was 
always 93.7 — 98.6% of the theoretical one) and above all 
due to the Ag electrode, which was applied to one or two 
faces and whose thickness and Young's modulus could 
not be controlled. The anelastic spectra are displayed in 
terms of the real part of the reciprocal Young's modulus 
or compliance, s — E^^ = s' — is", and of the elastic 
energy loss coefficient = s" / s' . Since the anelastic 
spectroscopy experiments were made in vacuum, some 



oxygen loss occurred when 800 K were exceeded, and the 
sample became darker. In most cases, however, we did 
not note any effect on the anelastic and even dielectric 
spectra; major changes occurred if 1000 K were exceeded 
with samples prepared without precautions against Bi 
loss during sintering. 



III. RESULTS AND DISCUSSION 
A. Undoped NBT 

1. The anelastic spectrum of NBT 

Figure [1] shows the compliance curves s' (T) and the 
losses Q"^ (T) of undoped NBT measured on heating and 
cooling at 1.5 kHz (thick black and grey lines); the thin 
black curves were measured at 8 kHz during heating and 
are practically coincident with those at 1.5 kHz, espe- 
cially in the real part. This means that the anomalies in 
s' are not due to relaxation of domain walls or defects, 
but they reflect the intrinsic elastic constants plus the fast 
relaxation of the order parameters involved in the struc- 
tural transformations, which can instantaneously follow 
the sample vibration. Some frequency dependence is ob- 
served in the losses, and may be due to the relaxation of 
domain walls. There are sharp peaks at Ti = 823 K and 
T2 — 564 K (during heating) in the complex compliance, 
whose sharpness is in contrast with the broad anomalies 
in the other physical properties, including the acoustic 
ones in the GHz range, and with the structural changes, 
which are faint and spread over wide temperature ranges. 
The anomalies appear at slightly lower temperature dur- 
ing cooling, with hystereses of 3 and 13 K for Ti and T2, 
respectively. 



2. Phase transitions in NBT 

The picture of the structural transformations in pure 
NBT is still confused and contradictory; therefore, in or- 
der to relate the elastic anomalies of Fig. [1] to struc- 
tural changes, we replotted our data in Fig. [2^), together 
with the dielectric susceptibility and the temperature de- 
pendencies of various structural parameters measured in 
other experiments. 

For clarity, the labels of the ordinates are omitted; only 
the temperature dependences are relevant here. Starting 
from the bottom panel a), the curves s' and mea- 
sured during heating determine the temperatures Ti and 
T2. In the next panel (ordinates: < e' < 3000) curve 
1 is the dielectric permittivity e' measured during heat- 
ing on poled NBT at 10 kHz. The step at - 470 K 
occurs at the same temperature where the macroscopic 
polarization vanishes, and the P — E hysteresis loops 
become narrow or pinched in the center, and which is 
therefore called depolarization temperature T^; it de- 
fines the transition between ferroelectric (FE) and possi- 
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FIG. 1: Reciprocal Young's modulus s' = 1/E' and elastic 
energy loss of NBT measured during heating (black) and 
cooling (gray) at 1.4 kHz (thick lines) and 7.4 kHz (thin line, 
only heating). 



bly antiferroelectric (AFE)^"— state. The antiferroelec- 
tric nature of the phase at T > Td is not unanimously 
acceptedii"— and is deduced, besides from the shape of 
the P — E hysteresis curves, from the broad maximum 
of e' around T„i ^ 630 K, which would signal the tran- 
sition from AFE to paraelectric (PE). The width of this 
maximum, a possible slight frequency dispersion, espe- 
cially at high frequencies^^ the fact that no structural 
transformation at Tm but instead several indications of 
phase coexistence on nanoscale existj^ii^ii^ induce several 
authors to liken this diffuse transition to a relaxor transi- 
tion not only in the solid solutions with BT or other fer- 
roelectrics, but also in pure NBT4iiii2r— The frequency 
dispersion is however definitely smaller than for tradi- 
tional relaxorsi^ii^ Due to limitations in the temperature 
range of our apparatus, we could neither reach the T-C 
transition nor the maximum of e' at T™, and therefore 
we report curve 2 from Ref. ^, with a reduced step at Td 
and an additional step, indicated by an arrow, connected 
with the R-T transition at a temperature higher than our 
T2. The position, amplitude and shape of such a step is 
highly variable from experiment to experiment while 
a dielectric anomaly has never been observed at Ti. 

Figure HJ:) presents the temperature dependence of 
the rotations of the octahedra determined by neutron 
diffraction^'^^ (ordinates: < w < 8°). In the case 
of NBT, the tetragonal structure arises from rotations 
of the octahedra about their c axes (a^a^c^ in Glazer's 
notatioii^^) , giving rise to superlattice peaks from the 
M point of the Brillouin zone, while the rhombohedral 
structure is due to rotations about all three principal di- 
rections of the cubic cell, with alternate sign along each 
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FIG. 2: Temperature depenndece of various properties of 
NBT from our measurements and from the literature, as ex- 
plained in the text. 



axis (a^a^a^ in Glazer's notation), giving rise to R-type 
superlattice peaks. The closed symbols (black — heating 
and gray — cooling) are the square root of the intensi- 
ties of these superlattice reflections^^ which should be 
proportional to the squares of the angles uj of rotations 
of type M and R.'^'^ They are scaled in order to overlap 
with UJ measured from Rietveld refinements^, (open sym- 
bols). The curves from the two experiments^^ii are in 
reasonable agreement with each other and represent the 
rotational order parameters of the two transitions. On 
cooling from the C phase, the M rotation (squares) has 
an onset at Ti ~ 820 ± 5 K, no hysteresis between heat- 
ing and cooling and a precursor tail extending a hundred 
kelvin above Ti, at least in the experiment of Vakhru- 
shev et aZ/^ This behavior is close to what expected for 
the condensation of a soft M phonon mode at the zone 
boundary, acting as order parameter of the T phase. The 
softening of such a mode has been deduced from the in- 
tensity of the quasielastic neutron scattering at the M 
point but not observed by other techniques, to our 
knowledge. On further cooling, the R rotations (trian- 
gles) start being observed below ^ 610 K and the anipli- 
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tude of the M rotations decreases. On the basis of these 
curves, the onset temperature for the R-T transforma- 
tion cannot be defined with an uncertainty smaher than 
40 K and there is a hysteresis of 60 K between heating 
and coohng. Again, a softening of the R phonon mode in 
the T phase has been deduced from the intensity of the 
quasielastic scattering at the R point. 

The cation shifts with respect to the octahedra de- 
duced from the neutron Rietveld refinements^ are plot- 
ted in Fig. [5Jl). In the ferroelectric R phase both the 
Na/Bi and Ti (not shown) shifts are parallel to the pseu- 
docubic [111] direction; in the T phase, they are parallel 
to the tetragonal axis but with opposite signs, so that 
the structure is only weakly polar, in agreement with the 
nearly antiferroelectric properties in this temperature re- 
gion. It is not clear whether the jump in cation shifts 
near 680 K may be put in relation with the broad dielec- 
tric maximum at T^. The refinements of the neutron 
scattering data of Jones and Thomas'^ indicate regions 
of coexistence of C and T phases below Ti and R and 
T phases around T2 (volume fractions in Fig. Un- 
fortunately, the large temperature intervals between suc- 
cessive measurements did not allow to clarify the contro- 
versy on the R-T coexistence regions, as discussed by the 
same authorsi^ A later diffuse neutron scattering studyi^ 
indicates a non trivial transformation between R and T 
structures over an extremely broad range of temperatures 
with a residual modulation of the R phase along the c 
axis of the parent structure. A modulated R structure 
is found also in recent TEM experiments, ^'^^ which are 
interpreted in terms of a low temperature i?3c structure, 
which above 473 K starts twinning parallel to the [010] 
pseudocubic planes. On further heating, the density of 
Pnma twin planes increases and the modulation of the 
polarization direction gives rise to antiferroelectric prop- 
erties until at 563 K the modulation disappears and the 
structure becomes homogeneous, identified as Pnma or- 
thorhombic rather than tetragonal. The temperatures 
of appearance and disappearance of the striation in the 
TEM images are indicated as vertical bars in Fig. 
and correspond to Td and T2. 



3. Identification of the elastic anomalies 

Based on Fig. [5] and the above discussion we can iden- 
tify the temperatures Ti and T2 of the peaks in the elastic 
compliance with the onsets of the T-C and R-T transfor- 
mations. The first case is obvious, in view of the coin- 
cidence of Ti with the onset of the M-type rotations of 
the octahedra (Fig. [2t) causing the tetragonal distortion. 
The fact that the T-C transformation is evident in the 
elastic but absent in the dielectric susceptibilityS^S con- 
firms that it is triggered by the non-polar rotational mode 
of the octahedra, as already discussed in connection with 
optical and high-frequency dielectric measurementSi^ 
The small hysteresis of 3 K between heating and cooling 
indicates that the transition is almost 2nd order. Accord- 
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FIG. 3: Left ordinates, continuous lines: dielectric permittiv- 
ity e' and losses tan 5 of NBT measured at 1 kHz on heat- 
ing (thick) and cooling (thin). Right ordinates, circles: elas- 
tic compliance s' and energy loss of NBT measured at 
1.5 kHz on heating (closed) and cooling (open). 



ing to the Landau theory of phase transitions without 
fluctuations, the expected effect on the elastic constants 
involved in the resulting tetragonal distortion is a nega- 
tive step when passing from cubic to tetragonal at Ti,— 
or equivalently a positive step in the compliances. The 
elastic anomaly at Ti has a consistent peaked component 
besides the step, but clearly corresponds to the almost 
2nd order T-C transition. The occurrence of the addi- 
tional peaked softening is common in perovskites and is 
usually explained in terms of fluctuation effects;^^ it will 
not be discussed further here. 

The R-T transition is characterized by a variety of ap- 
parently contradictory experimental results, which gen- 
erally agree on the difl'useness of the transition and broad 
range of coexistence of the two phases; the competition 
between the two M and R rotational instabilities even 
suggests the possibility of glassy behavior In some di- 
electric measurements there is a step near the region of R- 
T coexistence, but only on heating and not on cooling^i^ 
and with varying temperatures and shapeSi ^'^^i^^ The 
dielectric permittivity e' and losses tan 5 of one of our 
samples are plotted in Fig. [3] as thick (heating) and 
thin (cooling) lines, together with the compliance s' and 
elastic losses as closed (heating) and open (cooling) 
symbols. In spite of the limited maximum temperature, 
it is possible to see a step in the e' curves that corre- 
lates very well with the peak in s' . These compliance 
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curves s' and particularly the acoustic losses are sharply 
peaked at a temperature T2 with only 13 K of hystere- 
sis between heating and cooling and allow us to define 
a precise temperature for the R-T transition (note the 
sharpness of the losses at higher frequency at T2 in Fig. 
[1]). In the light of the most recent experiments,— 1^^'^^ the 
temperature T2 can be identified with the onset of the 
appearance of the R structure, which however appears 
gradually as a modulation between the M-type in-phase 
tilts of the octahedra about the c axis in the T phase and 
the R-type antiphase tilts of the R phasefi^ or as finely 
twinned R phase. Again, the fact that the anomaly is 
sharp and pronounced in the elastic compliance and less 
defined in the dielectric susceptibility suggests that the 
primary order parameter is the R-type tilting of the oc- 
tahedra, which in turn induces a change in the direction 
and amplitude of the cation shifts. Conversely, there is 
no trace of the ferroelectric transition at Td in the elastic 
compliance. This may also be due to the fact that the fer- 
roelectric transition does not seem to involve a change in 
the structure, but only in the size and correlation length 
of the domains. 

The present measurements do not allow yet to unify 
the different pictures of the R phase obtained from vari- 
ous techniques into a coherent one, but indicate that, in 
spite of the various indications of extreme heterogeneity 
on the local scale, short correlation lengths and relaxor- 
like behavior, the R-type instability responsible for the R 
phase produces an elastic anomaly comparable to those 
of other crystalline perovskites with normal structural 
transitions occurring at a well defined temperature. 

The above data can be compared with cn (T) mea- 
sured on a NBT crystal with the ultrasonic pulse echo 
technique,^' where a shallow minimum is found near 
580 K, with 25 K of hysteresis between heating and 
cooling. The minimum has been attributed to relaxation 
of the polar regions rather that to the phase transition 
itself, and, although its temperature is close to T2, it is 
much shallower than the peak in the reciprocal Young's 
modulus s' (T). Considering that the phase transition 
at T2 involves the rhombohedral distortion, it is possi- 
ble that the sharper response of the Young's modulus is 
mainly due to the shear elastic constants of C44 symme- 
try. Two Brillouin scattering experiments, which probe 
the acoustic properties at much higher frequencies, yield 
still different results. In one case^ a dip in sound veloc- 
ity and a peak in damping are found near 700 K, which, 
instead of coinciding with one of the temperatures of the 
T-C and R-T transitions, is between the two. The obser- 
vation has been explained in terms of strong fluctuations 
of the two coupled order parameters for the two T and 
R phases, plus disorder in the Na/Bi sublattice and has 
been put in parallel to relaxor ferroelectrics. In a more 
recent experiment,— the intensity of the quasi-elastic light 
scattering and a distortion of the phonon peaks from a 
lorentzian shape are found to be peaked at 550 K, which 
is about T2. 



B. Samples quality 

It is true that different techniques may produce differ- 
ent pictures of a same system; for example local probes 
like the NMR or EXAFS spectroscopies indicate larger 
and more disordered cation displacements than diffrac- 
tion techniques.^'' Yet, the quality of samples is impor- 
tant for NBT, since different samples may give different 
results from the same technique, notably in the dielectric 
and ferroelectric properties. The main issues are cation 
ordering, stoichiometry and related defects providing ion- 
ized electronic charges. Regarding cation ordering, there 
are ferroelectric perovskites like Pb(Sco.5Tao.5)03 (PST), 
where the level of Sc/Ta ordering can be controlled and 
produces marked effects on the temperature and diffuse- 
ness of the ferroelectric transitioui^^ In NBT, in some 
case superlattice | (111) peaks have been observed by x- 
ray diffraction^^ or TEM^Si^ and attributed to an alter- 
nation of Bi and Na along the three directions, but usu- 
ally cation order is found to be very weak and localized^ 
or absent)^ so that a random distribution is assumed 
The level of cation ordering in NBT and its solid solutions 
with other perovskites is not controllable as in PST, and 
prolonged annealings in order to improve ordering might 
rather cause loss of Bi,'''^ so that a spread of the NBT 
properties from sample to sample is not surprising. It is 
therefore appropriate to check whether the sharp elastic 
anomalies of NBT in Fig. [l]are the result of a sample of 
exceptional quality. 

Figure 2] shows the anelastic spectra of three samples 
of NBT prepared in different ways. Samples 1 and 2 were 
calcined for 1 h at 800 °C (sample 1) or 700 °C (sample 
2) and sintered for 2 h at 1150 °C with surrounding 
NBT packs in order to avoid Bi loss, while sample 3 was 
prepared with the same temperature schedules of sample 
1 but without precautions against Bi loss. Sample 2 is 
the one of Fig. [TJ 

It is apparent that the Bi/Na stoichiometry is an is- 
sue, since in sample 3 only a reminiscence of the T-C 
transition at Ti remains with the sign of the step in the 
compliance even reversed, and also the transition at T2 is 
smeared; in addition, strong irreversible relaxation phe- 
nomena occur above 750 K, which are absent in the sto- 
ichiometric samples. Such relaxation processes appear 
as an elastic energy loss two orders of magnitude larger 
than in the cubic phase of stoichiometric samples, and are 
presumably due to the migration of the light Na ion by 
a vacancy mechanism in the cation deficient Na/Bi sub- 
lattice. The effect is absent in a perfectly stoichiometric 
lattice. The dielectric susceptibilities of two samples pre- 
pared as sample 1 and 3 (not shown here) demonstrate 
that the less stoichiometric samples of type 3 are indeed 
of extremely bad quality also from the dielectric point of 
view, with enormous losses due to mobile charges, which 
totally conceal the transition at T^. On the other hand, 
samples of type 1 present e (T) curves entirely compara- 
ble to those of " good" samples in the literature; therefore, 
the sharp elastic anomalies at T2 should not be ascribed 
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FIG. 4: Reciprocal Young's modulus s' and elastic energy loss 
measured on heating (thick lines) and cooling (thin lines) 
at ~ 1.4 kHz on three samples of NBT prepared as described 
in the text. 

to an exceptional quality of those samples, but are a gen- 
eral property of NBT. In addition, the elastic anomaly at 
T2 clearly persists, although broadened, also in less stoi- 
chiometric samples of type 3, where the other transitions 
are washed out. This demonstrates that the rotational 
instability leading to the R-T transformation is also ro- 
bust against lattice disorder. 



C. Anelastic spectra of NBT-BT 

Figure [S] presents compliance s' and elastic energy 
losses curves measured on heating of a series of sam- 
ples with Ba substitution up to 6%. We do not include 
the curves for 8% for clarity and because we consider 
those measurements preliminary, due to the sample qual- 
ity. 

Upon substitution of Ba, the T-C transition lowers its 
temperature Ti (continuous arrows), broadens and fades 
away, until at a; = 0.06 it is no more discernible from the 
broad maximum that develops around Tme ~ 450 K. The 
R-T transition also lowers its temperature T2 (dashed 
arrows), and the corresponding elastic anomaly looses 
intensity and sharpness, but remains always clearly vis- 




T{K) 

FIG. 5: Reciprocal Young's modulus s' and elastic energy loss 
Q"^ of a series of (NBT)i_:r(BT)^ samples with < 2; < 0.06 
measured at ~ 1.4 kHz. The solid arrows indicate Ti, the 
dashed arrows T2 and the dotted T3 (Color online). 



ible. The peak in both the real and imaginary parts of 
the compliance maintains a rather sharp component at 
a decreasing T2, but its weight decreases in favor of a 
broad maximum at T^g which is predominant already at 
X = 0.03. The existence of two distinct anomalies at T2 
and True can be related to fact that there are at least three 
interacting but distinct order parameters, namely M-type 
and R-type rotations and cation shifts. Accordingly, the 
compliance presents three main anomalies: the one at Ti 
is obviously related to the M rotations responsible for the 
tetragonal distortion, while the other two should involve 
R-type rotations and cation displacements. The identifi- 
cation of the step with the R/T transition, and therefore 
the R-type rotations, is clear at a; = and x = 0.06, 
where the temperature T2 so defined coincides with that 
determined for the R/T transition by other techniques, 
and should therefore hold also for the intermediate tem- 
peratures (see also the phase diagram below). We re- 
main with the cation shifts as the main responsible for 
the maximum at T^e- 



1. The anomaly at T3 

At a; = 0.06 a new anomaly appears at a temperature 
indicated as T3 (dotted arrow) in Fig. [5] It would be 
natural to identify it with the ferroelectric transition at 
Td, but two observations cast doubts on this assignment. 
One is that deduced from the kink in e' and from 
the separation of the dielectric curves measured during 
heating and cooling is slightly different from T3. Fig- 
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FIG. 6: Left ordinates, continuous lines: dielectric permittiv- 
ity e' and losses tan5 of a sample of (NBT)i_i(BT)i with 
X — 0.06 measured at 1 kHz on heating (thick) and cooling 
(thin). Right ordinates, circles: elastic compliance s' and en- 
ergy loss measured at 1.4 kHz on heating; dashed lines: 
repetition of the measurement (heating and subsequent cool- 
ing) after one week of aging at room temperature. 



ure [6] compares the anelastic and dielectric curves of two 
samples of BNBT-6, where the correlation between the 
anomalies at T2 is excellent (hardly visible in s' but clear 
in tan S) but the identification of T3 with Td is dubious. 
The other observation is that in the elastic compliance 
of NBT there is no trace of Td (see Fig. [1]) and one 
would expect the same behavior at x = 0.06, unless the 
ferroelectric transition has changed nature. Indeed, at 
a; = the temperature Td separates a FE-R from a AFE- 
R phase, while at a; = 0.06 it separates the FE-T phase 
from a presumably AFE-C or pseudocubic one. Notice 
that the anomalies of both permittivity and compliance 
are evident during heating but disappear during cooling, 
indicating that in the highly disordered NBT-BT there 
is a slow structural evolution in the R phase, presumably 
domain coarsening. This is confirmed by the fact that 
the anomalies on heating require few days of aging at 
room temperature in order to fully develop. The dashed 
anelastic curves in Fig. [6] are the result of one week of 
aging at room temperature after a measurement at high 
temperature, and the two elastic anomalies are not yet 
fully developed. 
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FIG. 7: Real and imaginary parts of the dielectric permittiv- 
ity of a series of samples with < a; < 0.08 measured during 
heating at 10 kHz. 



D. Dielectric spectra 

The dielectric spectra have been extensively 
discussed,—"— and here we point out only the main 
conclusions and some differences with respect to the 
previous literature. 

In Fig. [7] a panoramic of the dielectric permittiv- 
ity and losses measured at 10 kHz during heating of 
(NBT)i_:r(BT)^ with < x < 0.08 is shown. The sam- 
ples with 2 and 3% Ba had rather high conductivities and 
hence losses, but the real part e' at 10 kHz, where the in- 
fluence of such conductivity effects is reduced, fit well the 
general trend of growing maximum of e' at and low- 
ering Td with increasing x. The step at T2 is invisible in 
this scale at a; = (but see Fig. [31) , it is indistinguishable 
from the ferroelectric transition due to the coincidence of 
T2 and Td for x > 0.02 but it is again separated from Td 
at a; = 0.06. A first difference between our measurements 
and those of Hiruma et al^ is that we see a single di- 
electric anomaly already at a; = 0.02 while they observe 
such a merging at x — 0.03. We also observed that tem- 
perature and width of these steps depend on the sample 
preparation and the polarization status, being generally 
sharper and at higher temperature in poled samples. 

Regarding the frequency dispersion of the dielectric 
maximum and the presence of a secondary maximum 
with even stronger frequency dispersion near Td^'^^ Fig. 
[5] shows the dielectric permittivity of an initially poled 
sample with x = 0.04 measured at different frequencies. 
There is some frequency dispersion at a; > 0.04, although 
of reduced magnitude with respect to what found in typ- 
ical relaxor ferroelectrics, but we never found the sec- 
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FIG. 8: Real and imaginary parts of the permittivity e of 
an initially poled sample of (NBT)i_^(BT)^ with x = 0.04 
(Color online). 
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FIG. 9: Phase diagram of (NBT)i_^(BT)^deduc6d from the 
present anelastic and dielectric (10 kHz) measurements. Grey 
lines: Tm, T2, Td from Ref. Isl . C = cubic, T — tetragonal, 
R — rhombohedral, F = ferroelectric, ~A — almost antifer- 
roelectric, P = paraelectric (Color online). 



ondary maximum near T^; we only find a neat step dur- 
ing heating in the poled state, with reduced sharpness 
and temperature during cooling, or heating in the un- 
poled state. The secondary maximum has been associ- 
ated with an incommensurate modulation between FE 
and AFE state j^^i^ as found in the highly ordered ferro- 
electric perovskite Pb(Ybo.5Tao. 5)0342 Incommensurate 
modulation of the R phase along the c axis of the parent 
T phase has indeed been found by neutron diffraction^^ 
and its persistence after the disappearance of the T phase 
has been suggested to be due to some type of ordering in 
the Na/Bi sublattice. Also TEM experiments indicate a 
R phase modulated by twins, ^'^^ but at T > Td. 

Our dielectric and anelastic spectra do not show fea- 
tures necessarily ascribable to incommensurate modu- 
lations; the anelastic losses remain high also below T2 
and Td, but this might also be due to generic domain 
wall motion. Apparently, such walls may give rise to 
incommensurate modulations, whose morphology likely 
depends on the sample quality. In this respect, we recall 
that the dielectric hump that is presumed to accompany 
incommensurate modulations is enhanced by introducing 
cation vacancies through La doping^ or in solid solu- 
tion with Ko.5Bio.5TiO3i2.-i4 or BaTiOs.^^ It is therefore 
likely that the exact stoichiometry and degree of ordering 
of cations determine the detailed structural distortions, 
their correlation lengths and domain morphologies, and 
the scarce consistency between the results of various ex- 
periments is not only due to differences in the properties 
that are probed but also in the quality of the samples. 



E. Phase diagram 

Figure [9] shows the phase diagram of 
(NBT)i_a;(BT)a;deduced from the present anelastic 
and dielectric (10 kHz) measurements during heating, 
together with that available in the literature^ (grey lines 
Trm T2, Td). The latter is mainly based on dielectric 
measurements and therefore does not include the border 
between T and C phases. It is sometimes assumed 
that the paraelectric phase above T^ is cubic, but this 
cannot be the case at low enough x, since the neutron 
diffraction data show that NBT is tetragonal above T^ 
and becomes cubic only at Ti = 820 K. 

The anelastic spectra of Fig. [5] allow us to draw a 
Ti (x) line that drops rapidly toward T^, roughly at the 
composition of the MPB; similarly to the NBT-PbTiOs 
phase diagramii^ It is not clear whether Ti merges with 
the Tm (x) line at x > 0.05 or dissolves into an aver- 
age cubic phase with local tetragonal distortions. The 
three open circles at a; = 0.05 refer to three different 
samples where the anomaly at Ti is extremely shallow 
and broadened. The region of the phase diagram defined 
by a; > 0.05 and Td < T < Tm is of uncertain nature. 
In (NBT)i_a;(BT)2; co-doped with Zr, with a particu- 
larly evident relaxor-like frequency dispersion in e' near 
Td, such a region has been defined as a frustrated region 
with coexistence of "multiple domain states" .— The fad- 
ing away of the anelastic anomaly at Ti with increasing 
X suggests that for x > 0.05 and T > Td the phase is 
pseudocubic with local tetragonal distortions whose cor- 
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relation length must be very small. 



The temperature T™ of the broad maximum of e' (T) 
is somewhat lower than in Ref. [s^, while the temper- 
ature T,„e of the maximum of the compliance is up to 
60 K lower than T„j. A similar phenomenon occurs in 
the relaxor ferroelectrics Pb(Mgi/3Nb2/3)i_a;Tia;03 and 
(Pb/La)(Zr/Ti)03, where the dielectric maximum is at 
a temperature higher than the maximum of the elastic 
compliance but also of the susceptibility deduced from 
NMR relaxation'*^''*® (but, unlike dielectric and anelastic 
spectroscopies, NMR probes also excitations with g > 0, 
see Ref. l49f ) . This has been explained as the result of the 
different sensitivity of the three techniques to the polar 
on one hand and non-polar rotational degrees of freedom 
on the other hand, all exhibiting a relaxational dynam- 
ics which is faster for the rotational ones.— Indeed, also 
the dielectric maximum in NBT-BT acquires a disper- 
sion in frequency with Ba doping. In NBT-BT, however, 
the R-type rotational degrees of freedom undergo a phase 
transformation at T2; therefore, the elastic maximum at 
T,ne should involve other rotational degrees of freedom, 
possibly also M-type rotations, which are coupled to the 
polar modes and hence yield dielectric and elastic sus- 
ceptibility maxima at similar temperatures. 



The T2 {x) line separating the R from the T phases 
joins Td [x) very soon, at cc = 0.03 according to Hiruma 
et al^ and already at x = 0.02 after our measurements. 
There is some variability in the values of T2 determined 
from the dielectric curves, depending on the polarization 
state, but overall there is good agreement with T2 de- 
termined from the anelastic spectra on unpoled samples. 
At X > 0.06 the Td (x) line starts rising, while T2 (x) 
continues to decrease, determining the MPB. 



As noted in the previous paragraphs, the dielectric 
anomaly at T2 may be not very sharp, especially in un- 
poled samples, and the marked anomalies measured at 
0.02 < X < 0.06 may mainly reflect the ferroelectric tran- 
sition. The compliance, however, is almost unaffected by 
the ferroelectric transition, and always presents a sharp 
step at T2. This means that the rotational instability 
of R type has a clear onset at a well defined tempera- 
ture T2 also at compositions where (NBT)i_a;(BT)a; dis- 
plays properties typical of highly disordered or relaxor 
systems. Notably, the sharp step in s' coexist with the 
broad maximum at Tme and a T-C transition that is al- 
most completely faded away at x > 0.05 (see Figs. [5] 
and [5]). It turns out, therefore, that the polar modes 
are strongly affected by disorder, characterized by short 
correlation lengths and undergo a diffuse or relaxor tran- 
sition; on the other hand, the R-type modes undergo a 
phase transition which is well-behaved from the elastic 
point of view, presumably with a long correlation length. 



F. Morphotropic Phase Boundary and possible 
monoclinic phase 

The composition-temperature phase diagrams of vari- 
ous solid solutions between a perovskite which is rhom- 
bohedral at low temperature and one which is tetra- 
hedral, have an almost vertical line that separates the 
two phases and is called morphotropic phase boundary 
(MPB). Among the most studied ferroelectric systems of 
this type are PZT^ and PMN-PT,^ After decades of 
studies on PZT, it has been found that an intermediate 
monoclinic phase exists in a narrow region at the MPBi^ 
The identification of such a monoclinic phase, where the 
direction of the polarization may continuously vary be- 
tween those in the R and T phases, was possible from 
high-resolution synchrotron x-ray diffraction experiments 
and careful analysis of the neutron diffraction data. The 
possibility for the polarization to continuously rotate be- 
tween the T and R phases through the monoclinic one is 
still debated r^i^ but is generally considered as the key 
factor that yields an easy change of polarization under 
stress and therefore high electromechanical coupling at 
the MPB compositions)^ 

The solid solution NBT-BT at room temperature 
passes from R to T at the composition x ^ 0.06, which 
is generally identified with the MPB, and where the elec- 
tromechanical coupling is maximum. Actually, there are 
few studies aimed at defining the MPB line T2 {x)r^^ 
which, to our knowledge, is defined by only two points 
at X = 0.06 and 0.07. Other tetragonal perovskites 
that give rise to a MPB with NBT are PbTiOg^ and 
Ko.sBio.sTiOs, but in the latter the composition of KBT 
at the border between R and T at room temperature has 
been reported to range between 0.17 and O.S.*** 

Besides better defining the MPB line, it would be in- 
teresting to investigate whether also NBT-BT has an in- 
termediate monoclinic phase that allows continuous ro- 
tation of the polarization between the pseudocubic [100] 
and [111] directions. The identification of a monoclinic 
phase by diffraction would certainly be more difficult for 
(NBT)i_^(BT)^ than for PZT, due to the larger disorder 
caused by the coexistence of ions with different charges, 
Na+, Bi'^+ and Ba^"*" in the A sublattice with respect to 
Zr*+ and Ti'^^, which have the same charge. The signa- 
ture of the M phase should be more evident in the elastic 
compliance versus temperature: in fact, the reciprocal 
Young's modulus of PZT has a large peak in correspon- 
dence with the MPB, 57 that may be simply explained in 
terms of the usual Landau free energy expansion with a 
coupling term linear in strain and quadratic in the po- 
larization, if the polarization may rotate in a monoclinic 
phase instead of switching between R and T phases. 

The case of (NBT)i_j:(BT)a; is more complicated than 
PZT, since there is no single transformation between T 
and M, but two distinct processes: a R-T transforma- 
tion at T2 and a diffuse transformation around T^e and 
Tjn- The first is well behaved from the elastic point of 
view, with the expected sharp step in the compliance 
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at the R-type rotational instability.^ If monoclinic dis- 
tortions exist, they must be connected with the diffuse 
process. The fact that it is the main feature in the dielec- 
tric susceptibility indicates that it is due to polar modes, 
while the strong peaked response in the elastic compli- 
ance recalls the behavior of the monoclinic phase of PZT. 
At variance with PZT, however, there should not be a 
monoclinic phase, but only local monoclinic distortions 
with short correlation lengths, thereby producing a dif- 
fuse response. This is in agreement with recent NMRSI 
and diffuse x-ray scattering experiments?^ both reveal- 
ing additional displacements of the Na/Bi ions away from 
the rhombohedral (111) direction along (100), producing 
a local monoclinic symmetry. The existence of mono- 
clinic regions is also supported by a diffuse scattering 
experiment, "'^^ from which it can be deduced that they 
are shaped as platelets one or two cell thick and 20 nm 
wide. 



IV. CONCLUSIONS 

We measured the complex dielectric permittivity and 
elastic compliance of (NBT)i_:c(BT)a; with < x < 
0.08. The anelastic spectra contain anomalies at both 
the temperatures Ti and T2 of the tetragonal-cubic and 
tetragonal-rhombohedral transitions, plus a broad max- 
imum at a temperature T^e, which is up to 60 K lower 
then Tm of the the diffuse dielectric maximum. The 
Ti (x) line in the x — T phase diagram drops rapidly 
with X and can be followed up to x ~ 0.05, where the 
transition becomes so diffuse and weak to be hardly rec- 
ognizable. At higher Ba substitution, the high temper- 
ature structure should be pseudocubic with short range 
tetragonal distortions. Instead, the elastic anomaly at 
T2 is extremely sharp in pure NBT and remains distinct 
up to X = 0.06, demonstrating that there is a well de- 
fined onset of the R-type octahedral rotations producing 
the R structure, in spite of the numerous indications of 
extremely diffuse transition and structural heterogeneity. 
The dielectric permittivity in our samples does not have 
a secondary maximum with frequency dispersion around 
the ferroelectric transition at T^, as often reported; in- 
stead, we observe steps both at Td and T2 (which however 



coincide for 0.02 < a; < 0.05), whose exact temperature 
and sharpness depend on the sample state and polariza- 
tion. These features are common to samples prepared 
under different temperature and time schedules, if pre- 
cautions are taken in order to avoid Bi loss during sin- 
tering; otherwise, the anelastic and dielectric spectra in- 
dicate extremely blurred structural transformations and 
defective material, especially above 800 K, where Na dif- 
fusion becomes possible through the vacancy mechanism. 

The anelastic and dielectric anomalies should be con- 
nected to the non-polar rotational modes of the octahe- 
dra and polar cation displacement as follows. The M-type 
rotations giving rise to the T phase becomes unstable at 
Ti , but the correlation length of such modes rapidly de- 
creases with Ba substitution, so that at x > 0.05 no long 
range T order exists any more. The R-type rotations 
giving rise to the R distortion are unstable at T2, and 
maintain long range correlations up to a; = 0.06, which 
is the Ba composition of the morphotropic phase bound- 
ary. The diffuse dielectric maximum at Tm and elastic 
maximum at Tme should be due to polar modes, whose 
correlation length remains short unless an external pol- 
ing field is applied or aging of the order of several days 
at room temperature is allowed. 

With Ba substitution, the broad peak in the elastic 
compliance becomes preponderant over the R-T transi- 
tion and, by analogy with the behavior of PZT with mon- 
oclinic structure at the MPB, the possibility is discussed 
that on the local scale the polarization may have an in- 
termediate direction between R and T. 
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